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ABSTRACT 

Context. Reflected light from giant planets is polarized by scattering, offering the possibility of investigating atmospheric properties 
with polarimetry. Polarimetric measurements are available for the atmospheres of solar system planets, and instruments are being 
developed to detect and study the polarimetric properties of extrasolar planets. 

Aims. We investigate the intensity and polarization of reflected light from planets in a systematic way with a grid of model calcula- 
tions. Comparison of the results with existing and future observations can be used to constrain parameters of planetary atmospheres. 
Methods. We present Monte Carlo simulations for planets with Rayleigh scattering atmospheres. We discuss the disk-integrated po- 
larization for phase angles typical of extrasolar planet observations and for the limb polarization effect observable for solar system 
objects near opposition. The main parameters investigated are single scattering albedo, optical depth of the scattering layer, and albedo 
of an underlying Lambert surface for a homogeneous Rayleigh scattering atmosphere. We also investigate atmospheres with isotropic 
scattering and forward scattering aerosol particles, as well as models with two scattering layers. 

Results. The reflected intensity and polarization depend strongly on the phase angle, as well as on atmospheric properties, such as the 
presence of absorbers or aerosol particles, column density of Rayleigh scattering particles and cloud albedo. Most likely to be detected 
are planets that produce a strong polarization flux signal because of an optically thick Rayleigh scattering layer. Limb polarization 
depends on absorption in a different way than the polarization at large phase angles. It is especially sensitive to a vertical stratification 
of absorbers. From limb polarization measurements, one can set constraints on the polarization at large phase angles. 
Conclusions. The model grid provides a tool for extracting quantitative results from polarimetric measurements of planetary atmo- 
spheres, in particular on the scattering properties and stratification of particles in the highest atmosphere layers. Spectropolarimetry 
of solar system planets offers complementary information to spectroscopy and polarization flux colors can be used for a first charac- 
terization of exoplanet atmospheres. 
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1. Introduction 

Light reflected from planetary atmospheres is generally polar- 
ized. The reflection is the result of different types of scattering 
particles with characteristic polarization properties. Polarimetric 
observations therefore provide information on the atmospheric 
structure and on the nature of scattering particles that comple- 
ments other observations. Systematic model calculations are re- 
quired to interpret the available polarimetry from solar system 
planets and prepare for future polarimetric measurements of ex- 
trasolar planets. 

Scattering processes. Rayleigh scattering occurs on particles 
much smaller than the wavelength of the scattered light. This 
process produces 100% polarization for a single right angle 
scattering. Rayleigh scattering is much stronger for short wave- 
lengths because the cross section behaves like cr oc I /A 4 , and it 
favors forward and backward scattering, which are both equally 
strong. The blue sky in Earth's atmosphere is a well known ex- 
ample of Rayleigh scattering by molecules. 
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Aerosol haze particles with a size roughly comparable to the 
wavelength can produce strongly forward directed scatterings. 
Depending on the structure of the particle, a high (p > 90%) 
or low (p w 20%) fractional polarization results for a scattering 
angle of 90°. For example, the maximum polarization for scat- 
tering by optically thin zodiacal or cometary dust is not more 
than x 30% (e.g. Leinert et al. 1981; Levasseur-Regourd et 
al. 1996), while a polarization close to 100% is inferred for sin- 
gle scattering of haze particles in Saturn's moon Titan (Tomasko 
et al. 2008). 

Liquid droplets in clouds produce a polarization because of 
refraction and reflection, which can be particularly high (> 50%) 
for scattering angles of about 140° for spherical water droplets, 
corresponding to the primary rainbow (see e.g. Bailey 2007). 
Clouds made of ice crystals reflect and refract light in many dif- 
ferent ways, and no distinct polarization features like rainbows 
are expected, except locally, where ice crystals may have very 
similar structures. 

Multiple scatterings in planetary atmospheres randomize the 
polarization direction of the single scatterings and lower the ob- 
servable polarization significantly. Therefore the net polarization 
of the reflected light depends not only on the scattering angle and 
the properties of the scattering particles, but also on the atmo- 
spheric structure. For this reason it is not suprising that a large 
diversity of polarization properties exists for the solar system 
planets. 
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Observations. Venus shows a low (< 5%) negative polariza- 
tion, which is a polarization parallel to the scattering plane, for 
most phase angles . In the blue and UV, a rainbow feature with 
a positive polarization of several percent is present (e.g. Coffeen 
& Gehrels 1969, Dollfus & Coffeen 1970), indicating that the 
reflection occurs mainly from droplets in optically thick clouds 
(Hansen & Hovenier 1974). 

For the giant planets, only observations near opposition 
are possible with earth-bound observations. Near opposition 
the disk-integrated polarization is low because single back- 
scattering is unpolarized and multiple scattering polarization 
cancels for a symmetric planet. 

With disk-resolved observations of Jupiter, Lyot (1929) first 
detected that the Jovian poles show a strong limb polarization 
of order 5-10%. To understand this effect one has to consider a 
back-scattering situation at the limb of a sphere, where locally 
we have a configuration of grazing incidence and grazing emer- 
gence (for a plane parallel atmosphere) for the incoming and the 
back-scattered photons, respectively. Photons scattered upwards 
will mostly escape without a second scattering, and photons scat- 
tered down have a low probability of being reflected towards us 
after the second scattering, but a high probability of being ab- 
sorbed or undergoing multiple scatterings. Thus photons that are 
reflected towards us by two scatterings travel predominantely 
parallel to the surface. Because the polarization angle induced 
in a single dipole-type scattering process, like Rayleigh scatter- 
ing, is perpendicular to the propagation direction of the incoming 
photon, a polarization perpendicular to the limb is produced. 

Measurements at large phase angles (« 90°) for Jupiter with 
spacecrafts detected a polarization of « 50% for the poles while 
the polarization is much lower (< 10%) for the equatorial region 
(Smith & Tomasko 1984). The high polarization at the poles can 
be explained by reflection from a scattering aerosol haze layer, 
while the polarization at the equator is low because of reflection 
from clouds. Towards short wavelengths (blue) the polarization 
at the equator increases strongly, indicating that also Rayleigh 
scattering contributes to the resulting polarization. 

For Saturn the polarization is qualitatively similar to Jupiter 
with an enhanced polarization at the poles at short wavelengths 
(blue). In the red the polarization level of the poles is lower than 
for Jupiter (Tomasko & Doose 1984). 

Uranus and Neptune display a strong limb polarization along 
the entire limb (Schmid et al. 2006a; Joos & Schmid 2007). 
Albedo spectra (e.g. Baines & Bergstralh 1986) and the polar- 
ization indicate that Rayleigh scattering is predominant in these 
atmospheres. 

An interesting case is Saturn's moon Titan, which has a 
thick scattering layer of photochemical haze that produces a 
very high disk-integrated polarization of ~ 50 % in the B and 
R band (Tomasko & Smith 1982). More recently the Huygens 
probe measured the scattering and polarization properties of the 
aerosol particles in great detail during its descent through Titan's 
atmosphere (Tomasko et al. 2008). 

The observations show that Rayleigh scattering is an impor- 
tant polarigeneric process in atmospheres of solar system ob- 
jects, in particular for Uranus and Neptune, and for the equato- 
rial regions of Jupiter and Saturn. Besides Rayleigh scattering 
one has to consider the reflection from haze particles (aerosols). 
Scattering by small aerosol particles (d < A) may be approx- 
imated by Rayleigh scattering. For large particles, d > A, the 
strong forward scattering effect and the reduced polarization for 
right angle scattering cause significant differences when com- 
pared to Rayleigh scattering. 



Clouds dominate in the atmosphere of Venus, and at longer 
wavelengths (red) also in Saturn and Jupiter. The reflection from 
clouds produces only a low positive or even negative polarization 
signal in Venus, Saturn or Jupiter, typically at a level p < 5 %. 
In a first approximation one may therefore treat clouds like a 
diffusely scattering layer producing no polarization. 

Polarimetric measurements of stellar systems with known 
extrasolar planets were attempted, but up to now no convinc- 
ing detection of the polarized reflected light from an extrasolar 
planet has been made (Lucas et al. 2009, Wiktorowicz 2009). 
The deduced upper limits on the polarization flux from the close- 
in planet indicate that these objects are not covered with a well 
reflecting Rayleigh scattering layer. 

Model calculations. The classical theory for the analytic solu- 
tion of the multiple scattering problem is treated in the seminal 
work of Chandrasekhar (1950), from which the polarization of 
conservative (non-absorbing) Rayleigh scattering planets can be 
derived. Van de Hulst (1980) gives a comprehensive overview 
on theoretical work up to that time including many numerical 
model results. 

Schmid et al. (2006a) put together available model results 
useful for parameter studies of the polarization from Rayleigh 
scattering atmospheres. This includes the following model re- 
sults: 

- Phase curves for the disk-integrated intensity and polariza- 
tion for finite, conservative (no absorption) Rayleigh scatter- 
ing atmospheres for different optical thicknesses and ground 
albedos from Kattawar & Adams (1971), 

- the limb polarization at opposition for semi-infinite Rayleigh 
scattering atmospheres with different single scattering albe- 
dos derived from formulas and tabulated functions given 
in Abhyankar & Fymat (1970, 1971) and Chandrasekhar 
(1950), 

- the limb polarization at opposition for finite, conservative 
(no absorption) Rayleigh scattering atmospheres for differ- 
ent optical thicknesses and ground albedos from tabulations 
given in Coulson et al. (1960). 

For Venus detailed models for the reflection from clouds 
were developed, which demonstrate nicely the diagnostic po- 
tential of polarimetric measurements (e.g. Hansen & Hovenier 
1974). More recent modeling of the polarization from planets 
was performed mainly to analyze and reproduce polarimetric 
observations of Jupiter and Titan from spacecrafts (e.g. Smith 
& Tomasko 1984; Braak et al. 2002; Tomasko et al. 2008). 

Another line of investigation now concentrates on the ex- 
pected polarization of extrasolar planets. The Rayleigh and Mie 
scattering polarization of close-in planets was investigated by 
Seager et al. (2000). These calculations consider planets which 
are unresolved from their central star and the polarization signal 
is strongly diluted by the unpolarized stellar light. 

Stam et al. (2004) modeled the polarization of a Jupiter-like 
extrasolar planet with methane absorption bands for three spe- 
cial cases and presented polarization spectra and wavelength in- 
tegrated phase curves. Also monochromatic phase curves for a 
non-absorbing clear and a hazy atmosphere are available (Stam 
et al. 2006). Other studies determined the expected polarization 
from clouds of terrestrial planets (e.g. Bailey 2007) or the polar- 
ization of extrasolar analogs to Earth (Stam 2008). 

Despite all these models systematic model calculations are 
sparse in the literature. For finite Rayleigh scattering atmo- 
spheres, polarization phase curves have been calculated only for 
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few selected cases. No results are available for the limb polar- 
ization of atmospheres with finite thickness and absorption. 

It is the goal of this paper to present a grid of model re- 
sults for Rayleigh scattering models with absorption and to ex- 
plore the model parameter space in a systematic way. The results 
should allow a comparison with observations and provide a tool 
for their interpretation. Additionally effects of selected devia- 
tions from simple Rayleigh scattering models will be discussed. 

In the next section the paper describes our scattering model 
and the Monte Carlo simulations. Section 3 presents the results 
from a comprehensive Rayleigh scattering model grid covering 
the three atmosphere parameters: single scattering albedo a>, op- 
tical thickness of the Rayleigh scattering layer t sc , and albedo 
of the underlying reflecting surface As. In Sect. 4 we explore 
the effects of a mixure of isotropic and Rayleigh scattering, of 
particles with a forward scattering phase function, and of two 
polarizing layers. In Sect. 5 we discuss spectral dependences. 
Section 6 highlights some special cases and diagnostic diagrams 
which may be of particular interest for the interpretation of ob- 
servational data. A discussion and conclusions are given in the 
final section. Appendix A describes the tables with the numeri- 
cal results of our calculations of intensity and polarization phase 
curves for a grid of 333 model parameter combinations. These 
are available in electronic form at the CDS. 



Because of the limited spatial resolution of most observations it 
is very hard to exactly measure the polarization near the limb. It 
is much less difficult to evaluate a disk-integrated polarization or 
polarization flux and to estimate and correct the degradation of 
the observed value with respect to the intrinsic value with a sim- 
ulation of the observational resolution or point spread function. 
This approach is described in detail in Schmid et al. (2006a) for 
seeing limited polarimetry of Uranus and Neptune. 

Therefore we mainly discuss the intensity weighted polariza- 
tion (q r ) = (Q r )/I, which is the equivalent to the disk-integrated 
radial polarization J Q r {r)2n r dr normalized to the geometric 
albedo. The geometric albedo A g is the disk-integrated reflected 
intensity of a given model at opposition normalized to the re- 
flection of a white Lambertian disk. It corresponds to 7(0°) in 
our calculations. 

The radial polarization curves are qualitatively similar for 
most models. The shape of the intensity curve varies signif- 
icantly from limb darkening to limb brightening for different 
model parameteres and cannot solely be described by the geo- 
metric albedo. Additionally we choose the Minnaert law expo- 
nent k as fit parameter for the shape of the center-to-limb inten- 
sity curve. The Minnaert law for opposition is 7(r) = 7 r= o// 2 *~ 1 , 
where p.(r) = (1 -r 2 ) 1/2 . This yields the following one-parameter 
fit curve I(r) = 7 r=0 (l - r 1 )^ 1 . 



2. Model description 



2.2. Atmosphere parameters 



Our planet model consists of a spherical body of radius R, illu- 
minated by a parallel beam. This geometry is appropriate for not 
rapidely rotating planets with a large separation, d » R, from 
the parent star. Each surface element is approximated by a plane 
parallel atmosphere. This simplification is reasonable for planets 
without an extended, tenuous atmosphere. 

2. 1 . Intensity and polarization parameters 

The intensity and polarization of the reflected light is described 
by the Stokes vector I = (7, Q, U, V). The linear polarized inten- 
sity or polarization flux is defined by the parameters Q - Iq- ho 
and U = 745 - 7135, where the indices stand for the polarization 
direction with respect to a specified direction in the selected co- 
ordinate system. In this paper only processes producing linear 
polarization are studied and therefore the Stokes parameter V 
for the circular polarization is omitted. We express the fractional 
polarization by the symbols 



Q 
1 ' 



u 
7' 



p = 



Ve 2 + u 2 



and the polarized intensity 

P i= ve 2 + u 2 . 



(i) 



(2) 



For the study of the limb polarization in resolved solar sys- 
tem planets at opposition we introduce the radial Stokes param- 
eter Q r , which is positive for an orientation of the polarization 
parallel to the radius vector r (perpendicular to the limb) and 
negative for an orientation perpendicular to r. The Stokes U r pa- 
rameter is the polarization direction +45° to the radius vector 
(see e.g. Schmid et al. 2006a for an illustrative description of the 
radial polarization). The polarization fraction is represented by 
q r and u r . 

The radial polarization curves q r (r) and Q r (r) can only be 
observed if the planetary disk is well resolved. The measured ra- 
dial profile depends strongly on the achieved spatial resolution. 



The plane parallel atmosphere is assumed to consist of a homo- 
geneous scattering layer that is either semi-infinite or finite with 
a reflecting (cloud or ground) Lambertian surface layer with a 
surface albedo. The basic model atmospheres are described by 
three parameters: 

- the single scattering albedo w, 

- the (vertical) optical thickness for scattering, t sc , of the scat- 
tering layer, 

- the albedo As of the surface below the scattering layer. 

The single scattering albedo o> is defined by the ratio be- 
tween the scattering cross section cr and the sum of absorption 
cross section k and scattering cross section <x, with the cross sec- 
tions multiplied by the fractions of scattering or absorbing parti- 
cles (/ sc or / abs ) 



AbsK + fsc<T 



(3) 



The value u — 1 indicates pure scattering (no absorption) while 
cj = is the other extreme of no scattering and just absorption 
(e.g. black dust). Similarly, a surface albedo of As =0 corre- 
sponds to a black surface, while a perfectly white Lambertian 
surface is defined by As = 1 . 

The optical depth for scattering t sc follows from the column 
density Z of the scattering layer: t sc = Z • cr, where cr is the 
scattering cross section per particle. The semi-infinite case cor- 
responds to t sc = 00 • We treat absorption like an addition of 
absorption optical depth to a layer with a given scattering opti- 
cal depth t sc , which is equivalent to reducing the single scatter- 
ing albedo. This approach is suited for discussing the reflected 
intensity and polarization inside and outside of absorption fea- 
tures like CH4 or HaO-bands, where k differs dramatically while 
cr is essentially equal. Then the total optical thickness t of the 
layer including absorption k is given by 



T = (fscCT + /abs*) ' Z = — . 



(4) 
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Fig. 1. Model geometry. The dashed line represents the trajectory of a 
reflected photon. 

The basic model grid (Sect. 3) considers only Rayleigh 
scattering (cr = o- Ray , t sc = T Ray ) as scattering process and 
Lambert surfaces with an albedo As below the scattering layer. 
Extensions, such as including non-polarizing isotropic scattering 
where cr = cr Ray + cr iso , haze layers or more than one scattering 
layer are discussed in Sect. 4. 

2.3. Geometric parameters 

The geometric parameters describe the location of the consid- 
ered surface point P and the escape direction of the photons 
(Fig. 1). A global coordinate system describes the orientation 
of the planet with respect to the star. Its polar axis is the sur- 
face normal at the sub-stellar point S', and the location of each 
point P is described by polar angle 5 and azimuthal angle (not 
drawn). 6 is also the photon's angle of incidence at point P. The 
escape direction, i.e. the location of the observer, is given by a 
polar angle a and azimuthal angles (not drawn), a is equivalent 
to the phase angle defined by the three (central) points: star or 
sun S, planet 0, and observer E (Earth). 

For the description of the scattering processes, a local coordi- 
nate system is set up at point P for the plane parallel atmosphere 
with surface normal z perpendicular to the planet surface in P, 
polar angle ft and azimuthal angle <p. 

In general, each point P can have individual atmospheric 
properties. Then the model outputs, the Stokes vector compo- 
nents /, Q and U, depend each on seven parameters: 

1(5, 0, a,x, t sc (S, 0), co(6, 0), A s (6, 0)) . 

This description allows calculation of the reflected intensity 
and polarization of each surface point on the illuminated hemi- 
sphere viewed from any direction. Obviously this large param- 
eter space needs to be simplified for a first parameter analysis. 
If we adopt the same atmospheric structure everywhere on the 
planet, t sc , oj and As are no longer functions of 6 and and we 
obtain a rotationally symmetric model geometry with respect to 
the line S - 0, which is independent of the azimuthal angle x- 

For extrasolar planets it will not be possible to resolve the 
disk in the near future. For disk-integrated results we can elimi- 
nate the dependence of the reflected intensity and polarization on 
the surface point parameters 5 and ft. Because of the rotational 
symmetry of the geometric model, the intensity and the polariza- 
tion then only depend on the polar viewing angle or phase angle 



a. Moreover the orientation of the polarization signal is either 
parallel or perpendicular to the scattering plane (the plane S-0- 
E), which we call the Q polarization direction. Q is defined pos- 
itive for a polarization perpendicular to the plane S-O-E and neg- 
ative for a polarization parallel to this plane. The U -polarization 
is zero in this coordinate system for symmetry reasons. 

For the full disk the integrated intensity and polarization sig- 
nals from a planet depend on the following parameters: 

I(a, t sc , oj, A s ) , Q(a, t sc , oj, A s ) . 

For solar system planets at opposition we obtain a rotation- 
ally symmetric scattering geometry (viewing direction is iden- 
tical to the axis of symmetry of the geometric model). We then 
have a scattering model which depends only on 6 or the normal- 
ized projected radius r = sin 6, and which is independent of 0. 
The resulting polarization will be in the radial direction either 
parallel or perpendicular to the radius vector r and therefore our 
model output is the radial Stokes parameter Q r (cf. Sect. 2.1). 
The Stokes U r parameter again has to be zero for a spherically 
symmetric planet. 

For exact opposition the dependences of the scattering model 
results can be described by the following parameters: 

l(r, t sc , oj, A s ) , Q r (r, t sc , oj, A s ) . 

These are the center to limb intensity curve and the center to 
limb radial polarization curves which both depend only on the 
atmospheric parameters. 

2.4. Monte Carlo simulations 

For our simulations we used the Monte Carlo code described 
in Schmid (1992), which was slightly adapted for the case of 
light reflection from a planet. Basically the code calculates the 
random walk histories of many photons in the planet model at- 
mosphere until the photons have escaped or are destroyed by an 
absorption process. After a sufficiently large number have es- 
caped, the scattering intensity and polarization of the reflected 
light can be established for different lines of sight. In our calcu- 
lations we assume that despite multiple scatterings the escaping 
photons emerge at the same point where they penetrated into 
the planet. In each scattering process the photon undergoes a di- 
rection and polarization change calculated from the appropriate 
phase matrix. The linear polarization of the photons in the sim- 
ulations is defined by the orientation y of the electric vector for 
the photon's electromagnetic wave. In a given coordinate system 
we can then evaluate the contribution to the Stokes intensity for 
each photon in Q oc cos 2y and U oc sin 2y direction. 

The escaping photons have to be collected in discret direc- 
tion bins (in our models phase angles a = 2.5°, 7.5°, . . . with a 
finite bin width Aa = 5°) to evaluate 1(a) and Q(a). These are 
then a mean photon intensity and polarization for that bin. Aa 
should be small to resolve any structure in 1(a) and Q(a), but 
also sufficiently large to collect enough photons for results with 
small statistical errors. The aim of our simulations is to reach at 
least the expected precision of observational data. The rotational 
symmetry imposed on our models helps to increase the bin size 
for the phase curve interval a^, which behaves like a^ oc sin a. 
This means that we have to divide the photon count per bin by 
the factor 27r sin aAa. The intensity is obtained by normalizing 
with the reflectivity of a white Lambertian disk. For a given sim- 
ulation the relative statistical errors (photon shot noise) are par- 
ticularly good for a « 90°, much less favorable for a = 2.5° and 
very bad for a = 177.5° where only a few photons will be col- 
lected, because the irradiated hemisphere of the planet is almost 
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invisible for this phase angle. For the center-to-limb curves we 
bin uniformly in 6 — arcsin(r) with a bin size of A6 = 5°, which 
requires an additional normalization by 2n sin 6 cos 6A6. 

The number of photons per model was chosen such that the 
number of reflected photons in phase angle bins relevant for ob- 
servations (a « 30° - 120°) are about N ~ 2- 10 6 when integrated 
over the whole disk. This corresponds to an error in polarization 
Ap = ^2/N = 0.1%. For the radial curves the total number of 
photons was increased such that the same precision was reached 
in most radial bins. No photons emerge at the exact phase an- 
gle a = 0°. Therefore for the limb polarization calculations we 
count all photons that are in the bin 0° < a < 5°, even though 
the calculation for the radial polarization includes the assump- 
tion that a = 0°. The error induced by this measure is smaller 
than the statistical error. 

A general guideline for the Monte Carlo technique for ran- 
dom walk problems is given in Cashwell & Everett (1959) and 
many Monte Carlo simulations for the investigation of light 
scattering are described in the astronomical literature (see e.g. 
Witt 1991, Code & Whitney 1995, Wolf et al. 1999). In Schmid 
(1992) a detailed description on many aspects of the employed 
Monte Carlo code are given; e.g. the general scheme of the code, 
the required transformations between the involved coordinate 
systems (star - planet, planet - plane parallel atmosphere, atmo- 
sphere - photon), the determination of the free path length, the 
treatment of isotropic scattering and Rayleigh scattering accord- 
ing to the Rayleigh phase matrix, an assessment of statistical 
errors, and a comparison with analytical calculations. 

3. Model results for a homogeneous Rayleigh 
scattering atmosphere 

This section discusses the model grid results for simple homoge- 
neous Rayleigh scattering atmospheres described by parameters 
0), t sc and As (cf. sec. 2.2) We discuss phase curves (Sect. 3.1) 
and radial profiles (Sect. 3.2) for selected cases and explore the 
full parameter space for disk-integrated results at a = 90° and 
a = 0° (Sect. 3.3). 

Many of the general dependences of these model results 
on atmospheric parameters were already discussed in previous 
studies mentioned in the introduction (Sect. 1). Compared to 
these our calculations are much more comprehensive and the ex- 
tensive model grid results are provided in electronic form (see 
Appendix A). An overview of the dependence of observable 
quantities, such as intensity, fractional polarization and polarized 
intensity, on atmosphere parameters is presented in diagrams 
which may be useful for the interpretation of observational data. 

The results presented in this section are in very good agree- 
ment with the previous calculations in Kattawar & Adams 
(1971), Stam et al. (2006) and Schmid et al. (2006a). 



Table 1. Reflectivity 1(a), polarization fraction q(a) and polarized in- 
tensity Q(a) phase curves for a very deep (t = 30) conservative (a> = 1) 
Rayleigh scattering atmosphere above a perfectly reflecting Lambert 
surface (surface albedo A s = 1). This model approximates well a con- 
servative, semi-infinite Rayleigh scattering atmosphere. Additionally 
the fit parameter a(a) for the parametrization of the polarized intensity 
Q(a) (Eq. 5) is given for relevant phase angles. 



an 


1(a) 


q(a) [%] 


Q(a) 


a(a) 


2.5 


0.795 


0.0 


0.0000 




7.5 


0.785 


0.4 


0.0031 




12.5 


0.766 


1.1 


0.0084 




17.5 


0.740 


2.1 


0.0155 




22.5 


0.708 


3.4 


0.0241 


1.85 


27.5 


0.671 


5.1 


0.0342 


1.86 


32.5 


0.630 


6.9 


0.0435 


1.87 


37.5 


0.587 


9.1 


0.0534 


1.89 


42.5 


0.542 


11.4 


0.0618 


1.91 


47.5 


0.497 


13.9 


0.0691 


1.94 


52.5 


0.453 


16.6 


0.0752 


1.98 


57.5 


0.410 


19.3 


0.0791 


2.03 


62.5 


0.368 


22.0 


0.0810 


2.08 


67.5 


0.329 


24.6 


0.0809 


2.14 


72.5 


0.292 


27.0 


0.0788 


2.21 


77.5 


0.259 


29.1 


0.0754 


2.29 


82.5 


0.228 


30.7 


0.0700 


2.39 


87.5 


0.199 


31.9 


0.0635 


2.49 


92.5 


0.174 


32.5 


0.0566 


2.62 


97.5 


0.150 


32.5 


0.0488 


2.77 


102.5 


0.130 


31.8 


0.0413 


2.95 


107.5 


0.111 


30.5 


0.0339 


3.16 


112.5 


0.094 


28.6 


0.0269 


3.42 


117.5 


0.079 


26.2 


0.0207 


3.76 


122.5 


0.066 


23.4 


0.0154 




127.5 


0.054 


20.3 


0.0110 




132.5 


0.043 


17.0 


0.0073 




137.5 


0.033 


13.7 


0.0045 




142.5 


0.025 


10.4 


0.0026 




147.5 


0.018 


7.3 


0.0013 




152.5 


0.013 


4.4 


0.0006 




157.5 


0.008 


2.0 


0.0002 




162.5 


0.005 


0.0 


0.0000 




167.5 


0.002 


-1.4 


0.0000 




172.5 


0.001 


-1.9 


0.0000 




177.5 


0.000 









The statistical error of the Monte Carlo calculation for 1(a) is 
smaller than 0.001 for all a. The uncertainty of the polarization 
fraction is less than 0.1 % for phase angles between 5 and 165 de- 
grees. Extrapolating the intensity / towards a = 0° with a quadratic 
least-squares fit to the first four points (a = 2.5°, . . . , 17.5°) yields 
a value /(0°) = 0.7970. This agrees with the exact solution 7(0°) = 
0.7975 from Prather (1974) to the third digit. 



Hulst (1980), and Bhatia & Abhyankar (1982) published a po- 

3.1. Phase curves larization curve for Rayleigh scattering in graphical form, but no 

tabulated values could be found in the literature. 

For the investigation of extrasolar planets, the phase dependence In Qur Monte Carlo simulation we tteat the ^mi-infinite at- 

of the disk-integrated polarization is of interest We discuss the here as Tsc = 30 and A x which ields esse ntially the 

phase curve for selected model cases (Fig. 2, left): a semi-mfimte „ a§ an infinite , bm ayoids infinite scatteri of 

and a finite scattering layer with different absorption properties SQme hotQns Qur of ^ cafje ^ tabulated ^ ^ j 
of the scattering and surface layers. 

The semi-infinite, conservatively scattering layer is a good 

reference case for an illuminated sphere and is often used for Intensity: The intensity phase curves 1(a) have their maximum 

scattering atmospheres. All irradiated light is reflected after one at a = 0°, when the whole illuminated hemisphere is visible, 

or several scatterings and the spherical albedo is equal to 1. An and they decrease steadily to zero at a — 180°, where only the 

intensity phase curve for isotropic scattering is given in van de dark side of the planet is seen. The intensity 7(0°) is equiva- 
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Fig. 2. Left: Phase dependence of the intensity /, fractional polarization q and polarized intensity Q for Rayleigh scattering atmospheres. Right: 
Radial dependence of the intensity /, radial polarization q r and radial polarized intensity Q r at opposition. Line styles denote: Semi-infinite case 
t sc = co (solid) for single scattering albedos to = 1 (thick), 0.1 (thin) and finite atmosphere (t sc = 0.3) with oj =1 (dashed) and 0.6 (dash-dot) for 
surface albedos A s = 1 (thick) and (thin). Also shown is the intensity curve for conservative semi-infinite isotropic scattering (dotted). 



lent to the geometric albedo. It is 0.7975 for the semi-infinite 
Rayleigh scattering atmosphere (Prather 1974), higher than for 
the semi-infinite isotropic scattering model (/i so (0°) = 0.690, van 
de Hulst 1980) or a white Lambertian sphere (/Lam(0°) = 2/3), 
because the Rayleigh scattering phase matrix favors forward and 
backward scattering. On the other hand the Rayleigh scattering 
intensity curve is lower for the range a « 52° - 120°. 

Of course, the reflected intensity decreases with absorption 
(with lower single scattering albedo co) in the atmosphere and 
with the albedo As of the underlying surface layer. The effect 
of absorption in the scattering atmosphere is important for thick 
layers, while the albedo of the underlying surface is important if 
the optical depths of the scattering region above is small. A quan- 
titative description of these dependences is given in Dlugach & 
Yanovitskij (1974) and Sromovsky (2005b). 

Polarization fraction. The disk-integrated polarization fraction 
q(a) is always zero for phase angles a = 0° and a = 180° 
because of the imposed rotational symmetry. The polarization 
maximum is near the right-angle scattering configuration a « 
90°. 

The polarization for the semi-infinite, conservative Rayleigh 
scattering layer reaches a maximum of q = 32.6% for a = 95°. 
For reduced scattering albedo, e.g. due to absorption in a molec- 
ular band, q(a) increases (see e.g. van de Hulst 1980). This hap- 
pens because absorption strongly reduces the fraction of multi- 
ply scattered photons in the reflected light which have random- 
ized polarization directions. If the absorption is very strong then 
the reflected light consists essentially only of photons that made 
one single Rayleigh scattering. The polarization phase curve 
then approaches the Rayleigh scattering polarization phase func- 



tion p(a) = (1 - cos 2 00/(1 + cos 2 a) with a polarization close to 
100 % at a = 90°. 

For finite scattering atmospheres the polarization fraction q 
also depends on the albedo of the surface layer. In the models 
discussed in this section the polarization is only produced in the 
Rayleigh scattering layer, while reflection from the surface layer 
is unpolarized. Therefore the resulting polarization is low for a 
high surface albedo and high for a low surface albedo (see e.g. 
Kattawar & Adams 1971, Stam 2008). This reflects the relative 
contribution of the polarized light from the scattering layer with 
respect to the unpolarized light reflected from the surface under- 
neath. 

The peak of the polarization curve is shifted towards large 
phase angles (a « 110°) for models with thin scattering layers 
and high surface albedos, as previously described by Kattawar 
& Adams (1971). At large phase angles (a > 90°), when only a 
planet crescent is visible, the fraction of scattered photons hitting 
the planet initially under grazing incidence is relatively high. For 
a thin scattering layer grazing incidence helps to enhance the 
probability for a polarizing Rayleigh scattering. For this reason 
the polarized light from the Rayleigh scattering atmosphere is 
less diluted by unpolarized light reflected from the surface at 
large phase angles and the fractional polarization is higher. 

Polarized intensity. The polarized intensity Q(a), which is the 
product of polarization q and intensity /, is zero at a — 0° 
and 180°, while the maximum of the phase curve Q(a) is near 
a as 65°, depending slightly on the model parameters. The max- 
imum value for the polarized intensity, considering the entire 
parameters space, is Q mdx = 0.0812 for the semi-infinite, con- 
servative Rayleigh scattering atmosphere at a = 65°. It seems 
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unlikely that another type of scattering process and model atmo- 
sphere can produce a higher polarized intensity. 

The polarized intensity decreases with increasing absorption, 
because the drop in intensity is stronger than the increase in frac- 
tional polarization. The polarization flux is a rough measure for 
the number of reflected photons undergoing one single Rayleigh 
scattering. Second and higher order scatterings also add to the 
polarized intensity, but only at a much lower level. Adding ab- 
sorption can only reduce the number of such scatterings and 
therefore diminishes the polarized intensity. 

A very important property of the polarization flux Q is that it 
does not depend on the albedo of the surface layer As (assumed 
to produce no polarization) below the scattering region. 



3.2. Radial dependence for resolved planetary disks at 
opposition 

For the interpretation of the limb polarization of solar system ob- 
jects close to opposition, we discuss the radial or center-to-limb 
dependence of the intensity 1(f), the radial polarization q r {r) and 
the radial polarized intensity Q r (r) (Fig. 2, right) for the same 
model parameters as for the phase curves in Sect. 3.1. 



Intensity: The radial intensity curve 1(f) shows a pronounced 
limb darkening in the semi-infinite conservative case. For a 
strongly absorbing atmosphere, e.g. within an absorption band, 
the /(r)-curve becomes essentially flat. Thus for an absorbing 
(and homogeneous) semi-infinite atmosphere limb brightening 
cannot be produced. For comparison the center-to-limb inten- 
sity curve 1(f) for isotropic scattering and for a Lambert sphere 
(1(f) = l/n(l - r 2 ) 1 / 2 ) are also shown. 

For finite scattering atmospheres with an optically thin layer 
the center-to-limb intensity curve can show a limb brightening 
effect. Limb brightening occurs for a highly reflective scatter- 
ing layer (high cj) located above a dark surface (low As), e.g. a 
thin aerosol layer or a methane-poor layer above the methane- 
rich absorbing layer (e.g. Price 1978). Limb brightening is ob- 
served in solar system planets in deep absorption band (e.g. 
Karkoschka 2001; Sromovsky & Fry 2007). Limb brightening 
is investigated in more detail in Sect. 3.2.1. 



Radial polarization fraction: The radial polarization fraction 
q r (f) is always zero in the disk center because of the symmetry 
of the scattering situation. For all cases the polarization increases 
steadily towards the limb and reaches a maximum value close to 
the limb between r = 0.95 and 1.0. The polarization q r (r) is 
always positive, which means a radial polarization direction or 
limb polarization perpendicular to the limb. 

It is important to note that the limb polarization decreases 
with decreasing single scattering albedo cj (more absorption) in 
contrast to the situation at large phase angles. This indicates that 
the photons producing the limb polarization are more strongly 
reduced by absorption than the reflected "unpolarized" photons. 

The explanation is that singly scattered (i.e. backscattered) 
photons do not contribute to the limb polarization, while re- 
flected photons scattered twice or a few times are responsible 
for the largest part of the limb polarization. Absorption implies 
that a larger fraction of escaping photons are singly-scattered 
and therefore unpolarized at opposition. Note however that for 
the semi-infinite atmosphere, the maximum radial polarization 
is not reached in the conservative case. A slightly lower scat- 
tering albedo (to « 0.95) mostly reduces the amount of highest 




Fig. 3. Radial intensity curves for different model parameters normal- 
ized to the central disk intensity, showing examples of limb darkening 
and limb brightening. The solid line is the calculated model and the 
dotted line the best fit with the Minnaert law. (a), (b) and (c) are for 
conservatively (a) = 1) scattering layers above black surfaces (A s = 0) 
with scattering layer thickness t sc =0.1 (a), 1 (b) and oo (c). (d) is 
a thin (t sc = 0.2), highly absorbing co = 0.2 scattering layer above a 
white surface (A s = 0). 



order scatterings and thus the polarization fraction is somewhat 
enhanced when compared to the conservative case (see Schmid 
et al. 2006a, and Fig. 7 in Sect. 3.3). 

The fractional limb polarization q r (r) for finite scattering 
layers depends strongly on the albedo of the underlying surface 

qA r ) is high for low As and low for high As like for large 
phase angles. A low surface albedo decreases the photons with 
multiple scatterings in the plane perpendicular to the limb, which 
are polarized parallel to the limb, thus enhancing the polarization 
in perpendicular direction. Therefore the limb polarization of a 
bright layer over a dark one can be even higher than for a semi- 
infinite atmosphere. This is discussed in more detail in Sect. 3.3. 

Radial polarized intensity: The radial polarized intensity 
Q r (f) = q r (r) ■ 1(f) increases with r from zero in the disk cen- 
ter to a maximum at r > 0.9 and then drops at the very limb. 
For semi-infinite atmospheres, Q just decreases at all radii with 
decreasing single scattering albedo cj. 

For finite atmospheres, the limb polarization flux Q r (f) de- 
pends only slightly on the surface albedo. Decreasing As from 1 
to can increase Q r (f) at most ~ 0.002 for some models, while 
for most models Q r (r) is virtually constant. This is similar to the 
case for large phase angles. 

3.2.1 . Limb darkening and limb brightening vs. limb 
polarization 

For a surface with a low albedo As below a thin scattering layer 
the limb can be brighter than the disk center, an effect that is 
generally called limb brightening. In principle this effect should 
be called "a central disk darkening", because the low surface 
albedo As does not brighten the limb. It only absorbs more light 
in the center of the disk, where a higher fraction of photons reach 
the absorbing surface because of their perpendicular incidence 
when compared to the situation of grazing incidence at the limb. 
Despite this fact we will retain the term "limb brightening" and 
consider the limb brightness on a relative scale compared to the 
brightness of the disk center. 
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The limb darkening and limb brightening effect can be 
parametrized to a first order using the Minnaert law I(r) = 
/ r =o(l _ r 2 ) k ~ 1/2 . The Minnaert parameter k determines the shape 
of the curve, k — 1 corresponds to Lambert's law, k = 0.5 to a 
flat intensity distribution I(r) = Iq and k < 0.5 to limb brighten- 
ing. The k parameter was determined by fitting Minnaert's law 
to all modeled intensity profiles, fixing the intensity at the center 
and excluding the outermost point where the formula diverges 
for k < 0.5. With the exception of some cases mentioned below, 
most profiles can be fitted adequatly. 

In Fig. 3 different examples of limb darkening and limb 
brightening are shown along with the best fit of the Minnaert 
law. Intensities are normalized to the central disk intensity. There 
are two types of limb brightening curves. For very thin atmo- 
spheres the maximal brightness is measured at the very edge of 
the planet, while for a moderate optical depth the intensity raises 
slightly up to a certain radius (e.g. 0.9R p ) and then drops very 
close to the limb. This second case cannot be fitted with the one- 
parameter Minnaert law and is approximated here by a relatively 
flat curve k « 0.5. 

The limiting case of a conservative semi-infinite atmosphere 
yields a Minnaert parameter of k « 0.9. For to going towards 0, 
k tends to a flat intenstiy distribution k = 0.5. For finite atmo- 
spheres there is a strong dependence on the surface albedo As. 
For a strongly absorbing atmosphere over a bright surface (to 
low, As high) absorption is more likely towards the limb (k > 1), 
for a bright atmosphere over a dark surface the opposite is true 
(k < 0.5). In the latter case the central disk intensity is very low. 

Similar to limb brightening, the limb polarization is also 
enhanced for a bright scattering layer over a dark surface. 
However, there are fundamental differences between these two 
effects. Limb polarization arises only for a polarizing process 
like Rayleigh scattering, while limb brightening occurs also 
for non-polarizing scattering processes like isotropic scattering. 
Additionally limb brightening is the stronger the thinner the up- 
per bright layer, while limb polarization requires a sufficiently 
thick scattering layer above the dark surface. Finally, limb polar- 
ization can also occur for cases of limb darkening, e.g. the semi- 
infinite, conservative atmosphere. Therefore, the two effects pro- 
vide complementary diagnostics of the vertical structure of the 
atmosphere. 

This section explores the full parameter space for sim- 
ple Rayleigh scattering atmospheres. We explore the parameter 
space by varying one of the three parameters to, t sc and As (cf. 
sec. 2.2) while fixing the other two. We study the resulting inten- 
sity 7(90°), polarization fraction <?(90°), and polarized intensity 
g(90°) (Figs. 4 to 6). 

The shapes of the model phase curves for the intensity 
1(a), fractional polarization q(a), and polarized intensity Q(a) 
look very similar for different model parameters (see Fig. 2). 
Therefore it is reasonable for a model parameter study to select 
the results for the phase angle a = 90°, considering them as 
representative (qualitatively) for all phase angles. A phase angle 
a = 90° is ideal for extrasolar planets because all planets will 
pass through this configuration twice during an orbit, regardless 
of inclination. 

The same type of parameter study is presented for the limb 
polarization of planets at opposition (a = 0°). For this we deter- 
mine disk-integrated (averaged) quantities for the intensity and 
radial polarization from the model results (Figs. 7 to 9). The inte- 
grated intensity 7(0°) is equivalent to the geometric albedo, (q r ) 
is the intensity weighted average of the fractional polarization, 
and (Q r ) the integrated polarized intensity on the same scale as 



the integrated intensity. These quantities are determined as de- 
scribed in Sect. 2.1. 

Figures 4 and 7 show the dependence on the single scatter- 
ing albedo to. For a given scattering optical depth r sc a reduction 
in to is equivalent to an enhancement of the absorption k in the 
scattering layer. Strong differences in k(X) occur in planetary at- 
mospheres for molecular absorptions (e.g. due to CH4 or H2O) 
inside and outside the band while cr is essentially equal. 

In Figs. 5 and 8 the Rayleigh scattering optical depths from 
t sc = 10.0 to 0.01 are plotted. This illustrates quite well the pos- 
sible spectral dependence from short to long wavelengths (left to 
right) for a Rayleigh scattering atmosphere. Since the Rayleigh 
scattering cross sections is proportional to 1//1 4 , it is possible 
that a planet has t sc = 4 at 400 nm and r sc = 1/4 at 800 nm. 

3.3. Parameter study for quadrature phase a = 90° and 
opposition a = 0° 

The effect of the albedo As of the surface below the Rayleigh 
scattering layer is shown in Figs. 6 and 9. 

General results from the Figures 4 to 9 are: 

- lowering the Rayleigh scattering albedo to always results in 
a lower intensity /, and lower polarized intensity Q or Q r , 

- lowering the Rayleigh scattering albedo to results in a higher 
polarization q at large phase angles. Contrary to this the frac- 
tional limb polarization q r is reduced for lower to, 

- lowering the Rayleigh scattering optical depth t sc produces 
a strong reduction in the polarized intensity Q or Q r in the 
optically thin case r < 2 and causes essentially no change in 
Q or Q r in the optical thick case t > 2, 

- lowering the surface albedo As lowers the intensity / and 
enhances the fractional polarization q or q r , 

- changing the surface albedo As does not change the polar- 
ized flux Q and hardly Q r . 

The most important difference between the limb polarization 
(q r ) and the disk-averaged polarization <7(90°) is their opposite 
dependence on the Rayleigh scattering albedo co (see e.g. the 
middle panels of Figs. 6 and 9). This occurs because the limb 
polarization at opposition is mainly caused by photons undergo- 
ing two to about six scatterings rather than just one. 

Another difference is the influence of t sc on the fractional 
polarization: q drops with t sc for bright non-polarizing surfaces 
and increases for dark surfaces. It is more complicated at oppo- 
sition: the limb polarization is highest if the dark ground elim- 
inates photons that would otherwise scatter twice perpendicular 
to the limb, but the atmosphere is still thick enough to produce 
many photons that escape having scattered twice parallel to the 
limb. The maximum possible limb polarization (q r ) = 5.25% is 
reached for 77 = 0.8, As = and to - I. 

From the variation of t sc shown in Fig. 5 it can be seen that 
the polarized intensity 2(90°) saturates above r = t sc • to > 2. 
Therefore 2(90°) cannot probe deep atmospheric layers. For the 
intensity and fractional polarization, an absorbing ground under 
a conservatively scattering layer can be noticed even at r > 10. 

The polarized intensity Q(a) consists mostly of photons un- 
dergoing just one single Rayleigh scattering. Therefore, Q is not 
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Fig. 4. Intensity, polarization and polarized intensity at quadrature as 
function of single scattering albedo oj for optical depths t sc = oo (solid), 
0.6 (dashed), 0.1 (dash-dot) and surface albedos A s = 1 (left), 0.3 (mid- 
dle), (right). 
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Fig. 7. Geometric albedo, disk-integrated radial polarization and polar- 
ized intensity at opposition as function of single scattering albedo oj for 
optical depths r sc = oo (solid), 0.6 (dashed), 0. 1 (dash-dot) and surface 
albedos A s = 1 (left), 0.3 (middle), (right). 
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Fig. 5. Intensity, polarization and polarized intensity at quadrature as 
function of optical depth r sc for single scattering albedos oj = 1 (solid), 
0.8 (dashed), 0.4 (dash-dot) and surface albedos A s = 1 (left), 0.3 (mid- 
dle), (right). 
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Fig. 6. Intensity, polarization and polarized intensity at quadrature as 
function of surface albedo A s for optical depths r sc = oo (solid), 0.6 
(dashed), 0.1 (dash-dot) and single scattering albedos oj = oj R!ly = 1 
(left), 0.8 (middle), 0.4 (right). 



Fig. 8. Geometric albedo, disk-integrated radial polarization and polar- 
ized intensity at opposition as function of optical depth t sc for single 
scattering albedos oj = 1 (solid), 0.8 (dashed), 0.4 (dash-dot) and sur- 
face albedos A s = 1 (left), 0.3 (middle), (right). 



<U = 1 



a)=0.8 



co=0.4 



0.8 

0.6 

g. 0.4 

0.2 
0.0 

£ 4.0 

A 

^ 2.0 
cr 
V 

0.0 
a 0.02 



0.01 



0.00 



* \ ^ 
* \ 


s 

; — \ 

--* ^ 


— — — * --, . 


■> 




/ 

y 






- -~ " t 




— — — — y 













1.0 0.8 0.6 0.4 0.2 1.0 0.8 0.6 0.4 0.2 

As Ag 



1.0 0.8 0.6 0.4 0.2 0.0 

As 



Fig. 9. Geometric albedo, disk-integrated radial polarization and polar- 
ized intensity at opposition as function of surface albedo A s for optical 
depths r sc = oo (solid), 1 (dotted), 0.6 (dashed), 0. 1 (dash-dot) and sin- 
gle scattering albedos oj = oj Riiy = 1 (left), 0.8 (middle), 0.4 (right). 
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Table 2. Best fit parameter b(a, t sc ) for the parametrization of the po- 
larized intensity Q (Eq. 5). 





K30°,t sc ) 


fc(60°,r sc ) 


fr(90\r sc ) 


6(120°, r sc ) 


0.1 


0.38 


0.39 


0.39 


0.44 


0.2 


0.63 


0.65 


0.66 


0.75 


0.3 


0.77 


0.80 


0.82 


0.90 


0.5 


0.99 


1.02 


1.04 


1.10 


0.8 


1.23 


1.22 


1.26 


1.23 


1.0 


1.32 


1.33 


1.32 


1.25 


2.0 


1.52 


1.48 


1.42 


1.27 


10.0 


1.59 


1.59 


1.45 


1.28 



changed by processes which happen deep in the atmosphere or 
by diffuse scattering on the surface. Q is only reduced if the num- 
ber of single Rayleigh scatterings are reduced, e.g. because there 
is only a thin Rayleigh scattering layer or photons are efficiently 
absorbed high in the atmosphere. 

We can approximate the polarized intensity Q by the follow- 
ing parametrization: 



Q(a, t sc , o>) = Q(a, oo, 1) • (1 - e - fl(ff)T -) • co h(a ^\ (5) 



where a(a) and b(a, t sc ) are fit parameters. Table 2 shows the 
best fit parameter b(a, r sc ), while Q(a, oo, 1) and a(a) are listed 
in Table 1. 

For optically thick Rayleigh scattering atmospheres the po- 
larized intensity Q depends only on the single scattering albedo 
co, and the parametrization reduces to: 

Q(a, co, Tsc > 2) * Q(a, oo, 1) ■ co h(a ^> 2> . (6) 
At quadrature this is 

2(90°, co, t sc > 2) * 0.060 • lo 1A5 (±0.002) . (7) 

For the limb polarization flux {Q r (0°)} the dependence on 
co is much steeper because both / and q r drop with decreasing 
co, as can be seen from the bottom panels of Fig. 7. For thick 
Rayleigh scattering layers, r > 2, the w-dependence of the limb 
polarization flux is 

2r(T sc > 2) * 0.022 • co 423 (±0.001) . (8) 



4. Models beyond a Rayleigh scattering layer with a 
Lambert surface 

The three parameter model grid discussed in Sect. 3.3 provides 
an overview on basic dependences of simple Rayleigh scatter- 
ing models. In this section we describe a few results for particle 
scattering properties different from pure Rayleigh scattering, or 
models with more than one polarizing scattering layer. 



for Rayleigh scattering planetary atmospheres 
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Fig. 10. Polarization of atmospheres with Rayleigh and isotropic scat- 
tering at 90° as a function of isotropic single scattering albedo a>j, nor- 
malized to the case of pure Rayleigh scattering to, = 1 or to, = 0. 
Plotted models are: semi-infinite atmosphere (thick solid), r sc = 0.3, 
A s = 1 (dashed), r sc = 0.3, A s = (dash-dotted), and t sc = 0.05 (thin 
solid). All models are without absorption, i.e. w,- + co r = 1. 

4.1. Atmospheres with Rayleigh and isotropic scattering 

Pure Rayleigh scattering is a simplification for planetary atmo- 
spheres. Already for Rayleigh scattering by molecular hydrogen 
one needs to account for a weak depolarization effect, because 
the diatomic molecule is non-spherical. Another depolarization 
effect for scattered radiation occurs in dense gas because colli- 
sions with other particles take place frequently during the scat- 
tering process. In addition aerosols and dust particles can also be 
efficient scatterers in planetary atmospheres and the net scatter- 
ing phase matrix differs from Rayleigh scattering and should be 
evaluated, e.g. by using the more general Mie theory. 

A simple way for taking such effects into account in a first 
approximation is to use a linear combination of the Rayleigh 
scattering and isotropic scattering phase matrices 

S = w-R+(l-w)-I, (9) 

where w = cr Ray /cr and 1 - w = cr iso /cr are the relative contri- 
butions of the Rayleigh scattering and isotropic scattering to the 
total scattering cross section cr = cr Ray + cr iso . Note that the single 
scattering albedo co and the scattering optical depth t sc now in- 
clude both the Rayleigh and the isotropic scattering cross section 
(cf. Sect. 2.2). 

Isotropic scattering is non-polarizing. If the scattering in the 
atmosphere is composed of both isotropic and Rayleigh scatter- 
ing, then the fractional polarization and the polarized intensity 
are reduced by isotropic scattering, while the intensity is compa- 
rable (cf. Fig. 2). 

Figure 10 shows the fractional polarization p(90°) at quadra- 
ture as a function of cr iso /cr for a few representative cases. In the 
single scattering limit the decrease is linear, the strongest devia- 
tion from a linear law is found for the semi-infinite atmosphere 
because of the large amount of multiple scatterings. A similar 
behavior is found for other phase angles, as well as for the radial 
limb polarization at opposition. 
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4.2. Forward-scattering phase functions 

The high polarization of Jupiter's poles and the disk-integrated 
polarization of Titan (e.g. Tomasko & Smith 1982; Smith & 
Tomasko 1984) has been explained by the presence of a thick 
layer of polarizing haze particles. The derived single scattering 
properties indicate strong forward scattering and Rayleigh-like 
linear polarization with maximal polarization close to 100% at 
about 90° scattering angle. Particles that satisfy this behavior are 
thought to be aggregates that are non-spherical and with a pro- 
jected area smaller than optical wavelengths (e.g. West 1991). 
We investigate the polarization properties of a planet with such 
a haze layer. The particle scattering properties are implemented 
as described in Braak et al. (2002) using a simple parametrized 
scattering matrix of the form 



F(#) = 



F n m F 12 m 

F 12(d) F n (&) 
F 33 (#) 
F 44 (&) ) 



(10) 




30 60 90 120 

Scattering angle # [°] 



150 



180 



Fig. 11. Probability density function p(&) for Rayleigh scattering 
(solid), Henyey-Greenstein function with asymmetry parameter g = 0.6 
(dashed) and g = 0.9 (dotted). 



where # is the scattering angle and 



Fn(&) = P H c(g,'&) = 



1- 



(1 + g 2 -2gcos#)( 3 / 2 > 



Fn(0) _ cos 2 § - 1 
fn(0) ~ Pb cos 2 0+1 



Fnift) 2cosi 



Fu(0) cos 2 #+r 



F 44 = 0. 



(11) 

(12) 

(13) 
(14) 



Fu(§) or Phg(&) is the Henyey-Greenstein phase function with 
the asymmetry parameter g (see e.g. Van de Hulst 1980). g = 
corresponds to isotropic scattering, g = 1 to pure forward scat- 
tering, g < to enhanced backscattering. Since haze particles 
have been shown to be strongly forward scattering, we limit our 
discussion to the two cases g = 0.6 and g = 0.9. 

Figure 11 shows the probability density function p(#) for 
Phg(&) m comparison with Rayleigh scattering. The probability 
density function for the scattering angle 1? is the phase function 
Fn(#) weighted by sin(i?) and normalized such that the inte- 
gral over p(#) equals 1. From this function the probability of 
the scattering angle within a certain interval is calculated by in- 
tegrating p(§) over this interval. One can see that for the haze 
models small scatttering angles (foward scattering) are greatly 
enhanced in comparison to Rayleigh scattering, while the prob- 
ability for backscattering is much lower. 

F 12(d) I F describes the fractional polarization of the 
scattered radiation as a function of the scattering angle. For scat- 
tering on haze particles it can be similar to Rayleigh scattering 
scaled by a factor p m , the maximal single scattering polariza- 
tion at 90° scattering angle. For a first qualitative analysis we set 
p m = 1 which is an upper limit that may slightly overestimate the 



resulting polarization. The other matrix elements are identical to 
Rayleigh scattering. 

Figures 12 and 13 show the phase and radial dependences 
for the haze models similar to the Rayleigh scattering case in 
Sect. 3.1 and 3.2. 



Intensity: The phase curves of the haze models differ from the 
Rayleigh scattering models mainly at small phase angles. The 
geometric albedo is lower for the haze models because backscat- 
tering is strongly suppressed compared to Rayleigh scattering. 
This is already discussed by Dlugach & Yanovitskij (1974), who 
calculated albedos for semi-infinite hazy atmospheres. Our cal- 
culations result in slightly higher albedos because of the inclu- 
sion of polarization. At phase angles around 90° the intensities 
are very similar for all models for non-absorbing atmospheres. 
An absorber greatly reduces the albedo of a planet with enhanced 
forward scattering, because many photons penetrate deeply into 
the atmosphere after the first scattering and then have a high 
probability of being absorbed. The radial intensity curves mainly 
reflect the lower geometric albedo, while the shape of the curve 
is similar for all models. 



Polarization fraction: The angle of maximal polarization is gen- 
erally larger for the haze models than for Rayleigh scattering. In 
the semi-infinite conservative case it is « 110° for haze as op- 
posed to « 95° for Rayleigh scattering. The shift to larger an- 
gles is particularly enhanced for a finite haze layer over a bright 
Lambert surface. However the maximal polarization decreases 
with increasing g. For strong absorption, both in or below the 
scattering layer, the polarization phase curve tends toward the 
single scattering function like in the Rayleigh case. The frac- 
tional limb polarization of haze layers can be much higher than 
for Rayleigh scattering layers, with disk-integrated values reach- 
ing (q r ) « 11% and peak maxima q r (r) « 20%. This is under- 
standable because the singly scattered (backscattered) photons 
which are unpolarized are strongly reduced for foward scatter- 
ing particles. 
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Fig. 12. Phase dependence of the intensity /, fractional polarization q and polarized intensity Q for a haze layer. Left: Semi-infinite case r sc = oo for 
single scattering albedos oj = 1 (thick), 0.6 (thin) and g = 0.6 (solid), 0.9 (dashed). The dotted line is the Rayleigh scattering case for comparison. 
Right: Finite atmosphere t sc = 0.3 with oj = 1 for surface albedos A s = 1 (thick), (thin) and line styles as for the left plot. 
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Fig. 13. Radial dependence of the intensity /, radial polarization q r and radial polarized intensity Q r at opposition for a haze layer. Left: Semi- 
infinite case t sc = oo for single scattering albedos co = 1 (thick), 0.6 (thin) and g = 0.6 (solid), 0.9 (dashed). The dotted line is the Rayleigh 
scattering case for comparison. Right: Finite atmosphere t sc = 0.3 with a> = 1 for surface albedos A s = 1 (thick), (thin) and line styles as for the 
left plot. 
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Fig. 14. Fractional polarization q as function of t sc of the upper layer 
at quadrature for models with Rayleigh scattering (solid) or isotropic 
(dotted) upper layer, and semi-infinite Rayleigh scattering lower layer 
with oj = 0.6 (thick) or Lambertian lower layer with A = 0.2 (thin), 
which provides the same reflectivity. 



Polarized intensity: The polarized intensity Q(a) is signifi- 
cantly lower for forward scattering phase functions than for 
Rayleigh scattering in the phase angle range a = 30° - 90° and 
for the limb polarization effect at opposition. It drops strongly 
with increasing g or increasing absorption. Like for Rayleigh 
scattering the polarized intensity is independent of the surface 
albedo As. The phase curves Q(a) show a shift of the maximum 
towards larger phase angles when compared to Rayleigh scatter- 
ing, in particular for models with thin scattering layers. 



4.3. Models with two polarizing layers 

Up to now we have treated the region below the scattering layer 
simply as a Lambert surface with an albedo As, which produces 
no polarization. In this section we explore model results for 
two polarizing layers with different absorption properties, where 
the lower layer can be a semi-infinite Rayleigh scattering atmo- 
sphere as described in Sect. 3.1. 

We focus on the question at what depth of the upper scatter- 
ing layer t sc the polarization properties of the underlying layer 
are washed out by multiple scattering and are no longer recog- 
nizable in the reflected radiation. 

Figure 14 compares the fractional polarization <7(90°) for 
three cases as a function of scattering optical depth of the up- 
per layer, t SCjU : a non-absorbing Rayleigh scattering layer above 
a semi-infinite, low albedo Rayleigh scattering atmosphere, the 
same scattering layer above a low albedo Lambertian surface, 
and an isotropic, non-polarizing scattering layer above the semi- 
infinite, low albedo Rayleigh scattering atmosphere. 

The reflected polarization shows no dependence on the po- 
larization properties of the underlying surface for deep scattering 
layers with t sc > 2. There are too many multiple scatterings to 
preserve this type of information from deeper layers in the es- 
caping photons. An imprint from the polarization of the lower 
layer becomes visible for thin scattering layers with t sc < 2. 
Particularly well visible is the polarization dependence on t sc 
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Fig. 15. Wavelength dependences of the model parameters total optical 
depth t(A) = r sc + t C h 4 of the upper layer, single scattering albedo o>(A) 
or o> u (A) of the upper layer, surface albedo A s or single scattering albedo 
u>i(A) of the lower layer. Two cases are considered: A layer of Rayleigh 
scattering with CH4 absorption above a white (As = 1) Lambert surface 
(solid), Rayleigh scattering layer without absorption o) u = 1 above a 
deep clear atmosphere with methane absorption (dotted). 



for the case where a polarizing layer is located below an isotrop- 
ically scattering layer. The polarizing lower layer only becomes 
apparent for t; so < 1. 

The same is true for the polarized intensity, because the re- 
flected intensity only shows a very weak dependence on the 
phase function of the scattering layer. The effects are also very 
similar for the limb polarization at opposition. 



5. Wavelength dependence 

The wavelength dependence of the reflected intensity and po- 
larization of a model planet can be calculated using wavelength 
dependent parameters t sc (A), co(A), and A$(A) or co\(A) for single 
or double layer models respectively. These parameters must be 
derived from a model with a given column density of scattering 
particles and mixing ratios for Rayleigh scattering and absorbing 
particles. 

As an example we selected parameters which approximate 
very roughly an Uranus-like atmosphere (e.g. Trafton 1976) con- 
sidering only Rayleigh scattering by H2 and He and absorption 
by CH4. In our first example we look at a homogeneous scatter- 
ing layer with methane absorption above a reflecting cloud layer 
with a wavelength independent surface albedo As = 1. This is 
a strong simplification for Uranus because the methane mixing 
ratio is of order 100 lower in the stratosphere than in the tro- 
posphere (Sromovsky & Fry 2007). Nevertheless it is a useful 
example for discussing basic effects of the wavelength depen- 
dence. 
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In a second example we make a first approximation for a 
methane mixing ratio that is varying with height, by having 
an upper layer of finite thickness without methane and a lower 
semi-infinite layer that includes methane. 

The Rayleigh scattering cross section of molecular hydrogen 
is given by Dalgarno & Williams (1962) as 



O" Ray,H 2 (/i) - — 



.14 • 10 13 1.28 • 10 6 1.61 



A 6 A* ' 

where A is in A and o-R ay ,H 2 in cm 2 /molecule. 
The total Rayleigh scattering optical depth is 



(15) 



^Ray — Cray 



(n, - If 

«H 2 - 1) 



2 ' 



(16) 



where Z, is the column density and «, the index of refraction of 
the z'-th constituent 1 . We use the same wavelength dependence as 
for the H2 cross section for all constituents. Our upper scattering 
layer has a column density Z = Z, = 500 km-am. 2 For the at- 
mospheric composition we adopt particle fractions of 0.5% CH4 
in the single layer case, and a methane free upper layer with 1 % 
CH4 in the lower layer in the two layer case. In all layers the He 
fraction is 15% and the rest is H2. 

Because of the strong wavelength dependence of the 
Rayleigh scattering cross section, t sc (A) changes significantly 



1 «h 2 = 1.0001384, n H e = 1.000035, « C h 4 = 1.000441 

2 1 km-am = 2.687 • 10 24 molecules cirr 2 



from the UV to the near-IR (Fig. 15, top panel). Keeping u) and 
As fixed (no absorber) yields the intensity and polarization re- 
sults given in Figs. 5 and 8 as function of t sc , which are in this 
case equivalent to results as function of A. 

The wavelength dependent single scattering albedo oj(A) fol- 
lows from the CH4 absorption optical depth tch 4 = Zch 4 kch 4 (^) 
and the Rayleigh scattering optical depth according to 



(0(A) = 



T SC W) 



1 sc 



(17) 



The absorption cross sections *cn,(^) were taken from 
Karkoschka (1994) and the resulting a)(A) is shown in Fig. 15. 

The intensity 1(A), fractional polarization q(A), and polar- 
ized intensity Q(A) is determined from the wavelength depen- 
dent model parameters for our two cases at quadrature and op- 
position (Fig. 16). 

At quadrature both examples show similar results. In both 
cases the polarization is enhanced and the polarized intensity is 
reduced within methane absorption bands, only the changes are 
less pronounced for a non-absorbing upper layer. The polarized 
intensity Q(A) also drops with wavelength, but it is overall higher 
in the second case because the polarizing Rayleigh scattering ex- 
tends to deeper layers. The biggest qualitative difference is seen 
in the continuum polarization q(A). In the case of an underlying 
reflecting cloud, q(A) drops towards longer wavelengths because 
of the smaller scattering optical depth above the diffusely scat- 
tering cloud. In the second case there is only polarizing Rayleigh 
scattering and no depolarization effect, so that the increasing ab- 
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sorption with wavelength in the lower layer results in a higher 
polarization. 

Similar spectropolarimetric models but with a Jupiter-like 
homogeneous atmosphere (higher column density, less methane 
than in our example) above both a dark surface As = and a re- 
flecting extended cloud were discussed by Stam et al. (2004) for 
a = 90°. The qualitative behavior of intensity and polarization 
with wavelength is quite similar to our example. However, for 
the same column density and methane fraction we find a signif- 
icantly lower intensity and higher polarization within methane 
bands. The origin of this discrepancy is unclear. A further com- 
parison with intensity calculations for a Neptune-like atmo- 
sphere by Sromosvky (2005a) shows a very good agreement at 
all wavelengths. Based on this we conclude that our model spec- 
tra should be correct. 

Intensity and polarized intensity at opposition behave qual- 
itatively similar to the large phase angle case. However the 
fractional polarization q r (X) is completely different. Absorbing 
particles in the upper layer tend to reduce the fractional limb 
polarization, while absorption in the lower layer enhances it. 
Observations of the limb polarization of Uranus and Neptune 
(Joos & Schmid 2007) show that the fractional polarization is in- 
deed enhanced within methane bands. Clearly for modeling limb 
polarization of these planets in absorption bands it is important 
to take into account the proper vertical stratification of the ab- 
sorbing component. A detailed model accounting for methane 
saturation and freeze-out to fit the observations is beyond the 
scope of this paper. 



6. Special cases and diagnostic diagrams 

We explore some special and extreme model cases in diagnostic 
diagrams of observational parameters for phase angle a = 90° 
and opposition. 



6.1. Fractional polarization versus intensity 

Figure 17 displays the diagnostic diagram for the reflectivity 
7(90°) and the relative polarization q(90°) at phase angle a = 
90°. Also indicated are the iso-contours for the polarization flux 
2(90°). 
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Fig. 17. Intensity vs. polarization at quadrature for the grid models. The 
shaded area indicates the possible range for Rayleigh and isotropic 
models. The symbols and lines indicate: semi-infinite conservative 
Rayleigh scattering (square), Lambert sphere (round), and black planet 
(diamond). The dash-dotted line shows semi-infinite models, the dashed 
and full line finite models without absorption (a> = 1) with surface 
albedo A s = 1 and respectively. The haze models shown in Fig. 12 
are indicated by crosses (high albedo) and plusses (low albedo). 
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Fig. 18. Geometric albedo vs. disk-integrated radial polarization at op- 
position for the grid models. The models are the same as in Fig. 17. 



The diagram shows points for special model cases and curves 
for the dependence on specific model parameters. The shaded 
area defines the area of observational parameters covered by 
our 3-parameter model grid for Rayleigh scattering (Sect. 3). 
Including isotropic scattering (Sect. 4.1) or having a vertically 
inhomogeneous atmosphere (Sect. 4.3) does not expand the cov- 
ered area. 

Figure 17 emphasizes that it is not possible to have a 
Rayleigh scattering planet with both very high albedo and polar- 
ization. A high albedo implies either a lot of multiple or isotropic 
scattering, which both reduce the fractional polarization. On the 
other hand a high polarization implies mainly single scattering 
and therefore strong absorption and a low reflectivity. The max- 
imal polarization at a fixed intensity is given by the model with 
a conservative (a> = 1) scattering layer over a dark (As = 0) 



surface and appropriate t sc . The semi-infinite atmosphere with 
varying co gives only slightly lower results than the former mod- 
els. The maximum of the product Q-ql- 0.060 is reached for 
the conservative semi-infinite atmosphere (t sc = oo). Since the 
polarized intensity Q is independent of the surface albedo As, a 
change in As is equivalent to a shift along the Q iso-contours in 
the diagram. 

The diagram also indicates the location of the haze models 
discussed in Sect. 4.2. Most of the haze models lie within the 
same area as Rayleigh scattering. Only for very thick haze layers 
with high single-scattering albedo is it theoretically possible to 
get somewhat higher fractional polarization for a given intensity. 

Figure 1 8 is the same diagram at opposition for the geometric 
albedo 7(0°), the disk-integrated limb polarization (q r ) and iso- 
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Fig. 19. Disk-integrated radial polarization at opposition vs. polariza- 
tion at quadrature for the grid. The models are the same as in Fig. 17. All 
Rayleigh scattering models lie in the dark shaded area, (partly) isotropic 
models also in the light shaded area. 
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Fig. 20. Disk-integrated radial polarized flux at opposition vs. polarized 
flux at quadrature for the grid models. The indicated models are the 
same as in Fig. 17. All Rayleigh scattering models lie in the dark shaded 
area, (partly) isotropic models also in the light shaded area. 



contours for the radial polarized intensity (Q r ). Like for large 
phase angles, the limb polarization for fixed intensity is highest 
for the conservative Rayleigh scattering layers over a dark sur- 
face. However, for t sc — > the polarization drops to 0%, while 
at large phase angles with As = it raises towards 100% when 
the few reflected photons are mainly singly scattered. The semi- 
infinite models provide a distinctly lower fractional polarization 
signal than a finite conservatively scattering atmosphere over a 
dark surface. The fractional limb polarization for very low albe- 
dos can be significantly higher for atmospheres with haze than 
for Rayleigh scattering, because the unpolarized backscattering 
is greatly reduced. 

For models with two polarizing layers with different absorp- 
tion properties, the results are located in the same area as for one 
layer above a surface. The limiting cases are models with a com- 
pletely dark lower layer (equivalent to a dark surface), and two 
identical layers (equivalent to a single semi-infinite layer). For 
atmospheres that contain also isotropically scattering particles 
the polarization is always lower and the intensity either slightly 
enhanced or reduced depending on a because of the different 
scattering phase functions. 



6.2. Polarization near a = 90° versus limb polarization 

For the prediction or the future interpretation of the polarization 
of extrasolar planets it is of interest to compare the polarization 
at phase angles near quadrature with the limb polarization at op- 
position. 

Figure 19 shows a diagram for the fractional polarization 
q(90°) and the fractional limb polarization (q r (0°)). Again the 
special models are indicated with black symbols and lines as in 
Figs. 17 and 18. 

A lower limit for the polarization fraction at phase angles 
a « 90° can be set from the limb polarization at opposition 
(q r ) for Rayleigh scattering or partly isotropic scattering atmo- 
spheres. For example a limb polarization of (q r ) « 2% implies 



a minimal polarization of q(90°) ~ 20%. The upper limit for 
the polarization fraction q(90°) is not well constrained by (q r ). 
The lower limit for Rayleigh or isotropic scattering may overes- 
timate the polarization at large phase angles only for very thick 
and bright haze layers. 

A tighter correlation is obtained for the polarization flux 
2(90°) and the limb polarization flux (Q r ), which is shown in 
Fig. 20. All Rayleigh scattering models are located in a nar- 
row area along a line from the origin (Lambert sphere / black 
planet) to the semi-infinite, conservative Rayleigh scattering 
model. Thus, for Rayleigh scattering atmospheres, one can pre- 
dict the large phase angle polarization flux Q(a) from the limb 
polarization flux (Q r ) and vice versa. The area is slightly broad- 
ened if isotropic scattering is included in the models, but the re- 
lation still holds quite well. Only very thick and high albedo haze 
layers show a significantly lower Q(90°) for a given (<2 r (0°)). 



6.3. Broadband polarized intensity 

Color indices of observational parameters are often relatively 
easy to measure and they are helpful for the characterization of 
atmospheres. From the atmosphere models they are obtained by 
averaging spectral results (Sect. 5) over the filter bandwidths. 

Here we discuss the colors for a Rayleigh scattering atmo- 
sphere with methane as a main absorber. It is investigated how 
the polarized intensity color changes as a function of methane 
mixing ratio and column density above a cloud or surface. Color 
indices are calculated by integrating Q(A) over the wavelength 
range foreseen for filters in the sphere/zimpol instrument (Beuzit 
et al. 2006). The filters are assumed to have flat transmission 
curves with cut offs at 555 and 700 nm (R-band) and 715 and 
865 nm (I-band). We concentrate on the color index of the polar- 
ized intensity Qi/Qr (Fig. 21). 

The polarized intensity is higher at shorter wavelengths, 
and Qi/Qr < 1 for all models because of the decrease of the 
Rayleigh scattering cross section with wavelength and the gen- 
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Fig. 21. Polarized intensity 2r as a function of color Qi/Qr where £>r 
is the broad-band R signal (555 to 700 nm) and Q\ the broad-band I 
signal (715 to 865 nm) for Rayleigh scattering planets with methane at 
90° phase angle. Indicated are models with constant methane mixing 
ratios / (dashed) or constant atmosphere column density Z in km-am 
(dotted) above a Lambert surface. 



eral increase of the absorption cross section of methane with 
wavelength. Qi/Qr is near 1 only for very thick atmospheres 
with very little methane or very thin atmospheres above a sur- 
face with wavelength independent scattering properties. In the 
former case Qr and Q\ are very high, in the latter very low. For 
intermediately thick atmospheres the color index Q\j Qr mainly 
depends on the methane mixing ratio, while Q R mainly depends 
on the column density. 

From this diagram we may predict that a color index of 
Qi/Qr ~ 0.25 - 0.5 could be typical for Rayleigh scattering 
atmospheres with methane absorption. Aerosol particles and ab- 
sorbers other than methane are expected to have a different spec- 
tral dependence of the scattering and absorption cross sections. 



7. Conclusions 

This paper presents a grid of model results for the intensity and 
polarization of Rayleigh scattering planetary atmospheres, cov- 
ering the model parameter space in a systematic way. The model 
parameters considered are the single scattering albedo a>, which 
describes absorption, the scattering optical depth of the layer t sc , 
and the albedo of a Lambert surface As. The results of these 
model calculations are available in electronic form at CDS (see 
Appendix A). In addition we explore models which combine 
Rayleigh and isotropic scattering, as well as particles with strong 
forward scattering and atmospheres with vertical stratification. 

Simple Rayleigh scattering models are a good first approx- 
imation to the polarization of light reflected from planetary at- 
mospheres because some amount of Rayleigh-like scattering by 
molecules or very small aerosol particles can be expected in 
any atmosphere. From the model grid, which basically provides 
monochromatic results, the spectropolarimetric signal can be 
calculated. This is done by considering the wavelength depen- 
dence of Rayleigh scattering and absorption in an atmosphere 
with given column density and particle abundance (see Sect. 5). 



The phase curves for the reflected intensity and polarization 
show a strong dependence on phase but they always have simi- 
lar shapes (see Fig. 2). However, the absolute level of the phase 
curve is a strong function of atmospheric parameters such as the 
abundance of absorbers or aerosol particles, the optical thickness 
of the Rayleigh scattering layer, or the albedo of the surface layer 
underneath (Sect. 3 and 4, see also Kattwar & Adams 1971). 

The model calculations demonstrate that polarimetric ob- 
servations would provide strong constraints on the atmospheric 
properties of the planetary atmospheres. An example is the po- 
larization flux Q(a) of the reflected light which for optically 
thick atmospheres is a simple function of the single scattering 
albedo roughly according to Q(a) oc cj b (b a 1.5). If both polar- 
ization and intensity can be measured, then one can distinguish 
between highly reflective and absorbing planets with or without 
substantial layers of Rayleigh-like scattering particles. 

According to the models a similar diagnostic is possible with 
observations of the geometric albedo, center-to-limb polariza- 
tion profile, and limb polarization for solar system planets near 
opposition. The limb polarization is in addition particularly sen- 
sitive to the vertical stratification of scattering or absorbing par- 
ticles located high in the atmosphere. 

The diagnostic potential is further enhanced if data for differ- 
ent spectral features, e.g. inside and outside of absorption bands, 
or from different spectral wavelength regions can be combined 
(Sect. 5, see also Stam et al. 2004). 

The calculations presented in this work are based on simple 
atmosphere models and they are therefore mainly useful for a 
first interpretation of data. For spectropolarimetric data of high 
quality, which are already available for solar system planets, one 
should use more sophisticated atmospheric models including a 
more detailed geometric structure, accurate abundances, and bet- 
ter scattering models for aerosol particles. With such models it 
might be possible for polarimetric studies to make a contribution 
to our knowledge on the rather well known atmospheres of solar 
system objects. 

Nonetheless the simple limb polarization models are of in- 
terest because they link the model results for large phase an- 
gles, suitable for extrasolar planet research, to models which can 
be easily compared with observations of solar system objects. 
Thus it may be possible to associate polarimetric observations 
of extrasolar planets to solar system objects. On the other hand 
it is possible to predict the expected polarization for quadrature 
phase of Uranus- and Neptune-like extrasolar planets with this 
simple model grid based on the existing limb polarization mea- 
surements of Uranus and Neptune (Fig. 20). 

Polarimetric measurements for extrasolar planets are ex- 
pected in the near future from high precision polarimeters. The 
measurements will first provide the polarimetric contrast, which 
is the ratio of the polarization flux from the planet Q(a) to the 
flux of the central star according to 



R 2 



R 2 



C(a) = — -q(a)-I(a) = — -Q(a), 



(18) 



where R is the radius of the planet and D the distance from its 
central star. 
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Very sensitive polarimetric measurements of stars with 
known close-in planets already exist, taken e.g. with the plan- 
etpol instrument (Hough et al. 2006; Lucas et al. 2009). This 
instrument measures the polarized intensity Q(a) of the planet 
diluted by the unpolarized flux of the central star. It is then dif- 
ficult to separate the fractional polarization q and the reflected 
intensity /. D is known from the radial velocity curve, but al- 
ready the radius of the planet R may be hard to derive if the sys- 
tem shows no transits. For photometrically very stable stars the 
reflected intensity 1(a) ■ R 2 /D 2 may be measurable with photom- 
etry of the phase curve with high precision instruments like most 
(e.g. Rowe et al. 2008). An uncertainty in the planet radius will 
affect the precision of the estimation of the normalized reflected 
intensity 1(a) (or reflectivity) of the planet. 

sphere, the future "VLT planet finder", which includes the 
high precision imaging polarimeter zimpol, could provide suc- 
cessful polarimetric detections (Beuzit et al. 2006; Schmid et 
al. 2006b). This instrument will be able to spatially resolve 
nearby (d < 10 pc) star-planet systems and allow a polarimetric 
search for faint companions to stars. In a first step only the dif- 
ferential polarization signal, i.e. the polarization flux Q(a), can 
be measured in the residual halo of the central star. The mea- 
surement of the intensity signal 1(a) might be possible if fur- 
ther progress in techniques like angular differential imaging is 
achieved. Even if a determination of 1(a) ■ R 2 /D 2 is possible, 
an uncertainty remains in the translation to normalized intensity 
1(a) if the radius of the planet is not known. 

Thus it may initially be quite difficult to measure intensity 
and radius. For this reason it is important to investigate the diag- 
nostic potential of the wavelength dependence in the polarization 
flux in more detail. For example the R-band and I-band yield a 
polarization color index Q\j £>r from which it should be possible 
to infer constraints on the Rayleigh scattering optical depth or 
the strength of absorption bands (see Sect. 6.3). Another route 
of investigation for atmospheres of extrasolar planets are mea- 
surements of the phase dependence of the polarization flux. For 
example the location of the maximum of Q along the phase curve 
is sensitive to the presence of aerosol particles, as discussed in 
Sect. 4.2. For planets in eccentric orbits the dependence of the 
polarization flux on the separation from the host star can be de- 
termined. 

One can hope that the current rapid progress in extrasolar 
planet observations continues, so that intensity measurements 
and accurate radius estimates for extrasolar planets become 
available soon after the first polarization flux detections, using 
the next generation of ground based telescopes and space instru- 
ments. Such instruments, if equipped with a polarimetric observ- 
ing mode, would allow a broad range of observational programs 
on the reflected intensity and polarization from planets. 
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Appendix A: Model grid tables 

Our extensive model grid of intensity and polarization phase 
curves for homogeneous Rayleigh scattering atmospheres 
(Sect. 3) is available in electronic form at CDS. Table A. 1 shows 
a sample of the first few lines and columns. The table is struc- 
tured as follows: Model parameters: Column 1: scattering optical 
thickness t sc , Column 2: single scattering albedo to, Column 3: 
surface albedo As, Model results: Column 4: spherical albedo 
A sp , Column 5: geometric albedo 7(0°), Column 6: limb polar- 
ization flux (Q r (0°)}, Column 7: 7(7.5°), Column 8: 7(12.5°),. . . , 
Column 40: 7(172.5°) Column 41: g(7.5°),..., Column 74: 
2(172.5°). Columns 7 to 74 are 7(a) and Q(a) spaced in 5 de- 
gree intervals. 7(0°) is equivalent to 7(2.5°) in our calculations. 
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Table A.l. Extract of model grid results. 





Ll) 


As 




1(0°) 


Q r (0°) 


7(7.5°) 


7(12.5°) . 


. 7(172.5°) 


2(7.5°) . 


• 2(172.5°) 


30.00 


1.00 


1.0 


1.0000 


0.7947 


0.02161 


0.7846 


0.7661 


. 0.0008 


0.00334 . 


. -0.00001 


99.00 


0.99 


1.0 


0.7947 


0.6378 


0.02108 


0.6290 


0.6130 


. 0.0007 


0.00331 . 


. -0.00001 


99.00 


0.95 


1.0 


0.5975 


0.4884 


0.01686 


0.4813 


0.4681 


. 0.0007 


0.00304 . 


. -0.00001 


99.00 


0.90 


1.0 


0.4794 


0.3980 


0.01316 


0.3918 


0.3807 


. 0.0006 


0.00278 . 


. -0.00001 


99.00 


0.80 


1.0 


0.3438 


0.2912 


0.00837 


0.2866 


0.2779 


. 0.0006 


0.00229 . 


. -0.00001 


99.00 


0.60 


1.0 


0.1966 


0.1707 


0.00341 


0.1676 


0.1623 


. 0.0004 


0.00146 . 


. -0.00000 


99.00 


0.40 


1.0 


0.1087 


0.0958 


0.00118 


0.0940 


0.0908 


. 0.0003 


0.00084 . 


. 0.00000 


99.00 


0.20 


1.0 


0.0470 


0.0418 


0.00024 


0.0410 


0.0396 


. 0.0001 


0.00039 . 


. 0.00000 


99.00 


0.10 


1.0 


0.0221 


0.0197 


0.00006 


0.0193 


0.0186 


. 0.0001 


0.00018 . 


. 0.00000 


10.00 


1.00 


1.0 


1.0000 


0.7949 


0.02141 


0.7848 


0.7662 


. 0.0008 


0.00338 . 


. -0.00001 


10.00 


1.00 


0.3 


0.8889 


0.7085 


0.02227 


0.6992 


0.6820 


. 0.0008 


0.00338 . 


. -0.00002 


10.00 


1.00 


0.0 


0.8833 


0.7042 


0.02236 


0.6950 


0.6779 


. 0.0008 


0.00336 . 


. -0.00002 


10.00 


0.99 


1.0 


0.8057 


0.6453 


0.02111 


0.6378 


0.6214 


. 0.0008 


0.00310 . 


. -0.00001 


10.00 


0.99 


0.3 


0.7875 


0.6326 


0.02112 


0.6238 


0.6076 


. 0.0007 


0.00342 . 


. -0.00003 


10.00 


0.99 


0.0 


0.7858 


0.6312 


0.02107 


0.6225 


0.6064 


. 0.0007 


0.00344 . 


. -0.00002 



2(2.5°), 7(177.5°) and g(177.5°) are very close to zero for all 
models and are not listed. 

All results are disk-integrated. Binning, normalization and 
errors are described in Sect. 2.4. For all calculations the number 
of photons was chosen such that A(Q/I) < 0.1% for phase angles 
a = 0° - 130°, and therefore (A7)/7 < 0.07% . 

The model grid spans the following parameters: t sc = 99, 
10, 5, 2, 1, 0.8, 0.6, 0.4, 0.3, 0.2, 0.1, 0.05, 0.01, co = 1, 0.99, 
0.95, 0.9, 0.8, 0.6, 0.4, 0.2, 0.1, A s = 1, 0.3, 0. Models for only 
three values of As are given because the polarized intensity is 
independent of As and the intensity drops nearly linearly with 
increasing As . Models with t sc = 99 were calculated only for As 
= 1 since the results are independent of As. Instead of a model 
with at = 1 and r sc = 99, we calculated the model with t sc = 30 
to reduce computation time, but the results are equivalent. 

The spherical albedo A sp in column 4 is the ratio of reflected 
photons in any direction to total incoming photons, while the ge- 
ometric albedo 7(0°) in column 5 is the disk-integrated reflected 
intensity at opposition normalized to the reflection of a white 
Lambertian disk. For our sample of Rayleigh scattering models, 
typically 7(0°) = (0.80 ± 0.06)A sp . 



